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LONG-TERM  GOALS 

The  long-term  goal  of  “Imaging  System  Performance  and  Visibility  as  Affected  by  the  Physical 
Environment”  is  to  investigate  the  relationships  between  physical  forcing  and  optical  properties, 
including  the  imaging  performance  parameters.  Our  efforts  are  aimed  at  supporting  the  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program,  which  is  devoted  to  the  topic  of  light  propagation  and  imaging 
across  the  air-sea  interface  and  within  the  surface  boundary  layer  of  natural  water  bodies. 

OBJECTIVES 

The  objectives  of  this  project  are  to: 

>  Determine  the  relationships  between  physical  processes  and  the  optical  properties  that  affect 
underwater  visibility  and  the  imaging  performance  parameters  (Chang  and  Twardowski,  2011). 

>  Investigate  the  effects  of  vertical  variability  (e.g.,  vertical  layering)  of  optical  properties  on  the 
image  performance  parameter,  the  modulation  transfer  function  (MTF). 


Here,  we  discuss  the  MTF  of  a  relatively  clear,  deep  water,  mesotrophic  environment  and  address  the 
following  issues: 

•  How  scattering  layers  and  their  relative  positions  from  detectors  and  targets  affect  the  MTF. 

•  To  what  degree  image  reconstruction  is  compromised  by  MTFs  computed  for  optically 
homogeneous  versus  optically  inhomogeneous  waters. 

APPROACH 

The  MTF  of  seawater  was  computed  using  profiled  measurements  of  inherent  optical  properties  (IOPs) 
conducted  between  the  near-surface  and  30  m  water  depth  in  the  Santa  Barbara  Channel  (SBC)  as  part 
of  the  ONR  RaDyO  project  (see  Chang  et  al.,  2010  and  Chang  and  Twardowski,  2011).  The  data  IOP 
processing  approach  utilized  factory-recommended  and  published  calibration  and  correction 
procedures  for  IOP  sensors  (e.g.,  Zaneveld  et  al.,  1994;  Pegau  et  al.,  1997;  Sullivan  et  al.,  2006;  Zhang 
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et  al.,  2009).  The  effects  of  stratification  and  scattering  layers  on  the  MTF  were  examined  for  two 
different  scenarios  (Figure  1): 

A.  A  virtual  imager  located  at  the  surface  and  virtual  targets  at  ranges  between  10  cm  and  30  m 
below  the  source. 

B.  Virtual  imagers  located  between  the  surface  and  30  m  and  the  targets  located  10  cm  below  the 
imagers. 


^n-t-2 
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Figure  1.  Schematic  diagram  of  two  scenarios:  A  (left).  Virtual  imager  (arrow)  and  virtual  targets 
(stars)  separated  by  scattering  layers,  z,„  over  a  total  distance,  Z,  and  B  (right).  Virtual  imagers 

located  a  distance,  zn,  above  virtual  targets. 


WORK  COMPLETED 

For  scenario  A  (Figure  1,  left),  the  collimated  decay  transfer  function  (DTF),  DC(VF),  and  MTF,  H(VP), 
at  range,  Z,  and  spatial  frequency,  \|/,  were  computed  following: 

Dc(y)s  =  c  -  S(\\i,  z)  dz 

where  the  s  denotes  the  scattering  layer,  c  is  the  beam  attenuation  coefficient  of  layer  s,  S  is  the  light 
scattered  back  to  the  acceptance  cone  of  layer  s,  and  an  integral  taken  between  the  imager  and 
distances  z/,  and  Z2  over  a  total  range  of  Z  to  the  target. 

0max 

S(v|/)  =  2 7i  Jo  J0  (2710 y)  P(6)  0  d0 

where  /3(d)  is  the  volume  scattering  function  (VSF),  dis  the  scattering  angle,  and  J0  is  the  zero-order 
Bessel  function.  The  near- forward  VSF  was  measured  using  a  LISST-100X  and  measurements  of  the 
beam  attenuation  coefficient  were  provided  by  an  ac-s.  The  focused  DTF  was  obtained  by 
transformation  in  terms  of  spatial  frequency  (27i0\|/)  prior  to  computation  of  the  MTF  (Wells,  1973): 

D(v)  =  J0  Dc(y  t)  dt 

H(v|/,z)s  =  e[Di(v)zi+D2(v)z2+-] 

Vertical  profiles  of  the  MTF  at  a  spatial  frequency  equal  to  1.0  rad'1  were  constructed  using  the 
equation  for  H(\|/,  z)s  above  by  summing  the  DTFs  between  the  surface  and  each  zn=  10  cm  layer  to  a 
maximum  range  of  Z  =  30  m  (Figure  2). 
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The  DTF  and  MTF  for  scenario  B  (Figure  1,  right)  were  calculated  following: 

D  Ok)  =  c  -  S(\| /) 

H(\|/,  z)  =  e'D(v)Z 

where  z  is  the  range  of  each  layer  (e.g.,  zn+2  to  Z),  with  similar  transformation  to  the  focused  DTF.  In 
case  B,  the  vertical  profile  of  the  MTF  for  \\i  equal  to  1.0  rad'1  was  constructed  with  multiple  virtual 
imagers  and  targets,  all  separated  by  the  same  distance  (10  cm)  (Figure  1). 

Scattering  layer  effects  on  the  MTF  for  the  two  different  scenarios  were  quantified  by: 

1)  Applying  a  power-law  fit,  y,  to  the  MTF  at  a  given  spatial  frequency  as  a  function  of  depth. 

2)  Determining  the  extinction  depth,  k,  of  MTF,  where  the  extinction  depth  was  defined  as  the 
depth  at  which  the  MTF  decreased  to  0.01  cycles/rad. 

RESULTS 

The  upper  water  column  MTF  at  small  spatial  frequencies  in  the  presence  of  scattering  layers  and 
stratification  can  be  orders  of  magnitude  less  than  the  MTF  in  mixed  conditions,  as  illustrated  in  Figure 
2.  The  distance  between  the  imager  and  target  contributed  greatly  to  the  magnitude  of  the  MTF  when 
imager  and  target  were  separated  by  a  scattering  layer. 


MTF|1*'  =  1)  MTF(*P  =  1) 


Figure  2.  Vertical  profiles  of  the  MTF(XF=  1  racT1)  shown  to  a  range  ofZ=  15  m  for  A:  Scenario  A 
and  B.  Scenario  B.  The  solid  dots  indicate  the  extinction  depth,  tc  Profile  color-coding  denotes 
specific  periods  of  the  SBC  RaDyO  experiment:  blues  =  strong  stratification,  yellow-green  =  upper 
water  column  mixing,  and  reds  =  strong  scattering  layers  from  localized  upwelling  (see  Chang  et  al., 
2010  and  Chang  and  Twardowski,  2011  for  more  details). 
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The  extinction  depth  determined  for  the  MTF(\|/  =  1  rad'1)  with  an  imager  at  the  surface  (scenario  A) 
did  not  vary  greatly  from  about  4.5  m  for  all  15  profiles;  the  minimum  value  of  k  was  2.9  m  in  the 
presence  of  strong  scattering  layers  and  maximum  k  was  4.8  m  during  mixed  conditions  (Figure  3). 
When  virtual  imagers  and  targets  were  placed  10  cm  apart  throughout  the  water  column,  k  varied 
between  4.5  m  (strong  scattering  layers)  and  9.8  m  (mixed  conditions). 

The  differences  in  y  between  scenarios  A  and  B  were  larger  during  stratified  conditions  (average 
difference  of  1.5)  as  compared  to  during  upper  water  column  mixed  conditions  (average  difference  of 
1.25).  The  smallest  difference  in  y  between  scenarios  A  and  B  (0.985)  occurred  during  periods  with 
strong  scattering  layers  (Figure  3). 


Prof  ile  Number 


Figure  3.  Upper:  MTF(XP=  1  rad1)  extinction  depth,  Middle:  Slope  of  the  power-law  fit  to  the 
vertical  profile  of  MTF(XF=  1  rad1),  and  Lower:  The  mixed  layer  depth  computed  as  the  depth  at 
which  the  cpg(660)  change  from  near-surface  cpg(660)  is  greater  than  0.05  in.  Yellow  =  scenario  A 
and  Purple  =  scenario  B.  The  profile  numbers  can  be  characterized  by:  1-7  =  strong  stratification, 

8-10  =  mixed,  and  11-15  =  strong  scattering  layers. 

These  results  imply  that  in  the  presence  of  near-surface  optical  scattering  layers,  an  underwater  image 
taken  from  just  below  the  surface  to  a  depth  of  a  mere  0.5  m  can  be  significantly  blurrier  than  if  the 
imager  was  located  tens  of  cm  below  the  surface.  Additionally,  analysis  or  reconstruction  of  an  image 
taken  by  a  near-surface  detector  of  a  target  located  at  just  a  couple  of  meters  below  the  surface  may  not 
be  possible  in  the  presence  of  strong  scattering  layers,  even  if  the  water  appears  relatively  clear  from 
the  surface. 

IMPACT/APPLICATIONS 

Knowledge  of  the  MTF  is  important  for  the  interpretation  of  images  from  underwater  electro-optical 
imaging  systems.  MTF-based  image  analysis  and  reconstruction  is  influenced  by  optical 
inhomogeneity  (i.e.  scattering  layers)  as  well  as  by  the  vertical  placement  of  the  target  versus  that  of 
the  detector,  or  imager  (Wells,  1973;  Frew  and  Voss,  1997).  Take,  for  example,  the  remark  made  by 
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Wells  (1973),  “You  can  lay  a  sheet  of  vellum  on  a  printed  page  and  read  through  it  quite  well,  but  the 
same  vellum  lifted  a  centimeter  off  the  page  completely  obscures  vision.”  This  statement  is  now 
outdated  as  vellum  is  no  longer  a  common  source  for  paper.  Here  is  another  way  to  consider  the 
problem:  one  should  not  feel  embarrassed  while  showering  several  feet  away  from  a  textured  shower 
door;  however  body  parts  could  be  fully  revealed  should  they  make  contact  with  the  door  (Figure  4). 

This  same  effect  applies  in  stratified  ocean  environments  where  a  target  and  imager  are  separated  by 
optically  distinct  layers.  These  layers  could  have  important  implications  on  imaging  performance  and 
analysis.  The  ability  to  determine  optimal  placement  of  passive  imagers  relative  to  target  locations 
would  constitute  dramatic  improvements  to  image  clarity  and  restoration  efforts.  Additionally,  image 
reconstruction  would  be  facilitated  with  knowledge  of  the  effects  of  scattering  layers  on  image 
distortion  and  imaging  performance  parameters. 


Figure  4.  These  photographs  of  Peso  the  dog  were  manipulated  in  Photoshop  to  simulate  Peso  well 
behind  a  textured  glass  shower  door  (Left)  and  Peso’s  nose  pressed  against  the  glass  (Right). 
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